Description

Implements the computation of the mesh characteristic properties using for_all loop over all GLL points and a scatterview that is required to store and access element-wise reductions. Implemented as a unified struct for both 2 and 3D, specfem::assembly::Info<dimension_tag> without specialization of the struct, but specialization of some of the methods under the hood.

Info Constructor Description:

The Info struct constructor unifies 2D and 3D mesh resolution analysis into a single templated definition. The constructor operates in two phases:

1. Element Processing Phase: Calls process_medium_elements<Dimension, Medium, Property>() for each medium/property combination in the assembly. This function uses a chunked domain iterator to traverse all GLL points, computing and reducing:
   - Material properties (vp, vs, rho) min/max
   - Domain bounding box coordinates
   - Element sizes (corner-to-corner distances)
   - GLL point spacing

Results are accumulated into an InfoScatters struct that holds all scatter reducers, including per-element views for velocity and distance bounds.
2. Derived Quantities Phase: A second parallel loop over elements computes the suggested timestep (from CFL condition using minimum GLL spacing and maximum velocity) and the largest minimum resolvable period (from average GLL spacing and minimum velocity).
These per-element values are reduced to obtain global bounds.

The InfoScatters struct uses dimension-templated ScatterMinMax views, allowing the bounding box reducer to be n-dimensional while scalar properties use 1D reductions.

Other necessary implementations

• Implements load_on_device for global_coordinates
• implements specfem::point::global_coordinates.coordinates() method which returns a view
• Moves mesh_entities::element.ngll member to public
• adds data_access type for global_coordinates
• adds vmin, vmax member functions for the different media adds tests for these
• adds point index copy constructor
• adds tests for scatterminmax struct
• adds tests for info computations (suggested dt etc.)

commit 1bcfcac4c3b88ff2a84e16715a216652bf905dfb
    Added 3D support

commit a0c2e1e4c632dfc318548bc3033160b2b5b97807
    Updated the code fully to be fully templated.

commit cc56d3d8edf5690d6ee08defb1ea982e5a68006c
    Updated the bounding box to use a registerarray

commit 535150ac4c6f9adad7e4dca927ebf57dc661bc14
    Updated the tests to match upstream fix

commit 3754efc1810817030e9f5566b7b3185916da548a
    Added ScatterMinMax struct for simpler interface of scatter reduction

commit 23c31ee6c9835b1d4d2a3ef69755b984dab01ba4
    Split out the core functions and tested them

commit 444cc8202a0766a55a0f7a519ea9eef59f92cd75
    Merge branch 'devel' into issue-1601

commit 94052e37ac2bfd6423ea65f21878627514dffd69
    Added all the missing files

commit 868eebb1b3af9e85a418d6aaa9e2f911a0a66a26
    Created 2D Info version.

commit 1132f0b89b117643e0b963fadee0e677424501d5
    Added computation of GLL point properties.

Issue Number

Closes #1601

Checklist

Please make sure to check developer documentation on specfem docs.

• I ran the code through pre-commit to check style
• THE DOCUMENTATION BUILDS WITHOUT WARNINGS/ERRORS
• I have added labels to the PR (see right hand side of the PR page)
• My code passes all the integration tests
• I have added sufficient unittests to test my changes
• I have added/updated documentation for the changes I am proposing
• I have updated CMakeLists to ensure my code builds
• My code builds across all platforms